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We investigated the Onsager relations in the context of electromagnetic 
constitutive relations of linear, homogeneous materials. We determined that 
application of the Onsager relations to the constitutive equations relating P 
and M to both E and B is in accord with Lorentz reciprocity as well as the 
Post constraint. 

1 Introduction 

In two seminal papers published in 1931 JH 13- with the assumption of micro- 
scopic reversibility, Onsager derived a set of reciprocity relations applicable to 
coupled linear phenomenons at macroscopic length scales. Fourteen years later, 
Casimir [3| improved the foundations of the Onsager relations. Initially consid- 
ered applicable to purely instantaneous phenomenons — or, at least, when "time- 
lag can be neglected" E p. 419] — the Onsager relations widened in scope as a 
result of the fluctuation-dissipation theorem [4| to time-harmonic phenomenons 
0. Sections 123-125 of the famous textbook of Landau and Lifshitz on statisti- 
cal physics provide a lucid introduction to the Onsager relations [6|, but we also 
recommend a perusal of a classic monograph by de Groot [7 1. A modern appraisal 
has been provided by Berdichevsky [8 1, whose paper motivated the work leading 
to this communication. 

Our focus is the correct application of the Onsager relations for linear electro- 
magnetic materials. This issue can be traced back to a 1973 paper by Rado 0. 
This paper contains a major conflict between a consequence of the assumption 



1 



of material response without any delay whatsoever and the Onsager relations as 
expounded by Callen et al. |5|. The former is definitely a noncausal assumption 
in electromagnetism ifTOl ITT1 . leading to false symmetries between the electro- 
magnetic constitutive parameters lfT2l . Furthermore, Rado considered E and H 
as primitive fields, but E and B are taken to be the primitive fields in modern 
electromagnetism H3J EH El- To the best of our knowledge, no other original 
investigation of the Onsager relations in electromagnetism exists. 

Due to the currently increasing emphasis on engineered nanomaterials lfT6lfT7l 
and complex electromagnetic materials ED ED, it is imperative that the applica- 
tion of fundamental principles (such as the Onsager relations) be carefully exam- 
ined with modern terminology. Accordingly, in the following sections, we first 
review the Onsager relations in general. Then we apply the Onsager relations to 
the electromagnetic constitutive relations of linear, homogeneous, bianisotropic 
materials. We show that a naive application to constitutive equations relating D 
and H to both E and B yields unphysical results, but that application to consti- 
tutive equations relating P and M to both E and B is in accord with Lorentz 
reciprocity [20] as well as the Post constraint ll2"Tll2"2l . 



2 Onsager relations 

Let us consider the linear macroscopic constitutive equations 

N 

L m = J2®mnF n , m £ [1, N] , (1) 

n=l 

where N > 1, L m are the Onsager fluxes and F m are the Onsager forces. The 
Onsager relations deal with the constitutive parameters $ mn . 

The derivation of the Onsager relations proceeds with the postulation of N 
state variables a n , n £ [1, N). The state variables are divided into two groups. 
The first N < N state variables are supposed to be even and the remaining N — N 
state variables are supposed to be odd with respect to a reversal of velocities of 
the microscopic particles constituting the linear medium; in other words, 

a m (t)a n (t + r) = a m (t)a n (t - r) , 

( m £ [1, N] andn £ [1, N] 
if <j or (2) 

{ m £ [N + 1, N] andn £ [N + 1, N] 
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and 

a m (t)a n (t + r) = - a m (t)a n (t - r) , 

( me [I, N] andn G [N + 1, AT] 
if or , (3) 

[ me[N+l, N]andn G [1, N] 

where the overbar indicates averaging over time t Q. 

In terms of the state variables, the Onsager fluxes are defined as 

d 

L m = —a m , m G [1, N] ; (4) 
the Onsager forces are defined as 

TV 

F m = - ^2 9m n a n , m G [1, N] ; (5) 

n=l 

and the coefficients g mn help define the deviation AS of the entropy from its 
equilibrium value as the quadratic expression 



I N N 



A5 = - EE^ 



^ m=l n=l 
I N N 



2 ^ t ^ ] QmnQ'rnQ'n • (6) 
m=7V+l n=iV+l 

In consequence of the microscopic reversibility indicated by © and ©, the 
constitutive parameters satisfy the Onsager relations 

<f> — $ 

f m G [1, N] andn G [1, iV] 
if < or (7) 

[ m G [iV + 1, N] and n G [iV + 1, N] 

and 

m G [1, N] andn E [N + 1, N] 

or . (8) 

m G [N + 1, iV] andn G [1, AT] 
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In an external magnetostatic field B dc , © and © are modified to 



and 



respectively. 



m G [1, iV] andn G [1, N] 
if { or (9) 

m G [# + 1, iV] andn G [iV + 1, N] 



m G [1, iV] andn e[N+l, N] 



if{ or , (10) 

m G [JV + 1, N] andn G [1, N] 



3 Application to Linear Electromagnetism 
3.1 Constitutive Equations for D and H 

Let us now consider a linear, homogeneous, bianisotropic medium. Its constitutive 
equations can be written in a cartesian coordinate system as 



Hj = J2l=i Cjk ° E k + Vj k o B k 



JG[1,3]. (11) 



We have adopted here the modern view of electromagnetism wherein E and B 
are the primitive fields while D and H are the induction fields lfT3lfT4lfT5ll . The 
operation o indicates a temporal convolution operation in the time domain, and 
simple multiplication in the frequency domain 

Now, D and E are even, but H and B are odd, with respect to time-reversal. 
With that in mind, we can rewrite (fTTT) compactly as 



N 



G[l, N), (12) 



n=l 



where F m = E m , F m+3 = B m , Q rn = D m and Q m+3 = H m for m G [1, 3]; 
furthermore, iV = 3 and N = 6. 
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With the assumption of microscopic reversibility, application of the Onsager 
relations © and ( fTUI) yields the following symmetries: 



A, 



'■mny-^dcj 

A mn (B dc ) 



Anm( B^cJ 



A 



B r 



nm\ "dc 
— A^^f — B 



dc) 



m e [1, 3] , n e [1, 3] 
m e [4, 6] , n G [4, 6] 
m e [1, 3] , n e [4, 6] 



Equations (TH5t imply that 



ejfe(Bd c ) 
^ifc(Bd c ) 
Cjfe(Bdc) 



-Cfcj(-B dc ) 



(13) 



(14) 



3.2 Constitutive Equations for P and M 

When considering a material medium, as distinct from matter-free space (i.e., 
vacuum), the presence of matter is indicated by the the polarization P = D — e Q E 
and the magnetization M = fi~ X B — H, where e D and [i are the permittivity and 
the permeability of matter-free space. Linear constitutive equations for P and M 
can be stated as 



Pj = ELi xfi ° E k + x ^ o B, 



(2) 



(4) 



J €[1,3], 



(15) 



where 



(i) 



l6 jk 
(2) 

_ (3) 

Xjk 



(i) 



(16) 



and Sjk is the Kronecker delta function. 

As P is even but M is odd with respect to time-reversal, we can rewrite (fl~5l) 

as 



N 



Rn 



*mn ^n, m G [1, JV] , 



(17) 



n=l 
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where R m = P m and R m +3 = M m for m G [1, 3]. As the microscopic processes 
underlying the constitutive parameters in (fT71) are reversible, \l/ mn must satisfy © 
and <H0l>; thus, 



^ 'mn (Brf c ) 



^ram( B^ c ) 



\1> fB 



del 



m G [1, 3] , n G [1, 3] ' 
m G [4, 6] , n G [4, 6] 
-^ nm (-B dc ) , m G [1, 3] , n G [4, 6] 



(18) 



whence the symmetries 



.(2) 



B 



<fcJ 



(3), 



del 



Xjk ( B <fc) = -Xjfcj (-B<fc 



(19) 



are predicted by the Onsager relations as the macroscopic consequences of micro- 
scopic reversibility. 



3.3 The Conflict 

Equations (fT9l imply that 



^'fe(Bd c ) = v k j(—~B dc ) 
Cifc(Bdc) = Cfcj(— Bd c ) 



(20) 



by virtue of (fl~6b . 

But d20lh disagrees completely with (O3. Let us reiterate that both (fT4lh and 
d^Ol c. come about from the application of the Onsager relations, contingent upon 
the assumption of microscopic reversibility. Yet, at most, only one of the two must 
be correct. 



3.4 Resolution of the Conflict 

Onsager's own papers help resolve the conflict. His papers were concerned with 
motion of microscopic particles, and he considered his work to hold true for heat 
conduction, gaseous diffusion and related transport problems. The Onsager forces 
must be causative agents, while the Onsager fluxes must be directly concerned 
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with particulate motion. This understanding is reinforced by subsequent com- 
mentaries (5J|7|. 

Therefore, in order to correctly exploit the Onsager relations in electromag- 
netics, we must isolate those parts of D and H which indicate the presence of a 
material, because microscopic processes cannot occur in matter-free space (i.e., 
vacuum). The matter-indicating parts of D and H are P and M. Hence, d^l ) 
must be accepted and d7?t must be discarded. 

With Brf c = 0, the symmetries (l20b coincide — unlike ([T4l — with those 
mandated by Lorentz reciprocity [1201 Eqs. 23]. Also unlike (TBI) , the symmetries 
(l20l are compatible with the Post constraint ll2Tll22l 

E & = E en (2D 

j=i i=i 

which must be satisfied by all (i.e., Lorentz-reciprocal as well as Lorentz-nonreciprocal) 
linear materials. These two well-known facts also support our decision to discard 
CD in favor of®. 

4 Concluding Remarks 

In this communication, we first reviewed the Onsager relations which delineate the 
macroscopic consequences of microscopic reversibility in linear materials. Then 
we applied the relations to the electromagnetic constitutive relations of homoge- 
neous bianisotropic materials. We determined that a naive application to consti- 
tutive equations relating D and H to both E and B yields unphysical results, but 
that application to constitutive equations relating P and M to both E and B is in 
accord with Lorentz reciprocity as well as the Post constraint. 
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